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£\j " Abstract 

A number of scalar invariant characterizations of the Kerr solution are presented. These 
, characterizations come in the form of quality factors defined in stationary space-times. A 

1 quality factor is a scalar quantity varying in the interval [0, 1] with the value 1 being attained 

if and only if the space-time is locally isometric to the Kerr solution. No knowledge of the 
Kerr solution is required to compute these quality factors. A number of different possibilities 
arise depending on whether the space-time is Ricci-flat and asymptotically fiat, just Ricci- 
fiat, or Ricci non-flat. In each situation a number of quality factors are constructed and 
analysed. The relevance of these quality factors is clear in any situation where one seeks a 
rigorous formulation of the statement that a space-time is "close" to the Kerr solution, such 
qq \ as: its non-linear stability problem, the asymptotic settlement of a radiating isolated system 

OO ■ undergoing gravitational collapse, or in the formulation of some uniqueness results. 
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1 Introduction 

The Kerr metric [16| is a fundamental exact solution in General Relativity providing the unique 
set of models describing the region outside a rotating black hole, so that its physical relevance is 
unquestionable. The uniqueness properties of the Kerr solution stem from powerful mathematical 
1 results stating, roughly speaking, that it is the only analytic, asymptotically flat vacuum solution 

of the Einstein field equations (without cosmological constant) representing a black hole with a 
regular event horizon (a precise, up-to-date formulation of this uniqueness result and its proof 
can be found in O [3]). The generalization of this result for non- analytic solutions of the field 
equations has been an open question for many years and recent progress towards its solution has 
been reported in [T4J, [TJ [15] 

A completely open problem regarding the Kerr solution is its stability under nonlinear perturb- 
tions. The non-linear stability of the Kerr solution — to be more precise, the non-linear stability 
of its domain of outer communication — can be formulated as the global existence of a solution of 
the vacuum Einstein equations corresponding to initial data whose maximal Cauchy development 
is close to a subset of the Kerr's domain of outer communication. Results in this direction already 
exist for other well known exact solutions to the Einstein field equations [8j [9] [TUJ [4] [TTl [20] and 
therefore, given the relevance of the Kerr solution, it is important attemping to obtain similar 
results for it. In this regard, it is crucial to establish what is the precise meaning of the statement 
that a solution of the field equations be "close to the Kerr solution" . 
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Another point of relevance of the Kerr solution arises from its status as the stationary equi- 
librium limit reached in the evolution of an isolated system undergoing gravitational collapse. 
Due to its uniqueness properties just mentioned, it is widely believed that after the gravitational 
wave content has been radiated away such systems will eventually settled down as the exterior 
region of a Kerr solution. Any rigorous proof of this statement requires to introduce the notion of 
"approaching the (exterior) of the Kerr solution" which in turn relies on a further notion of "being 
close to Kerr" . With the aid of this notion one could devise numerical simulations of a radiating 
isolated system and test whether in the asymptotic regime it is close to the Kerr solution. 

The purpose of the present work is to put forward a number of local and invariant criteria 
which enable us to answer the question of when a given space-time is "close" to the Kerr solution. 
By local and invariant we mean that the criteria are formulated in terms of scalar quantities, so 
that they have a perfectly defined value at each point of the space-time, and are independent of the 
coordinate systems used. This makes them suitable for their application to the problems discussed 
in the previous paragraphs. Our approach consists in introducing quality factors to measure the 
deviation of a given solution from the Kerr space-time. A quality factor is a dimensionless scalar 
quantity varying in the interval [0, 1] and such that the unity value is attained on a neighborhood 
of a point if and only if the space-time is isometric to the Kerr solution there. Ideally one seeks 
quality factors which are constructed from invariantly defined quantities of the geometry. The 
existence of such type of quality factors permits to quantify invariantly, by a normalized non- 
negative real number, "how close" to the Kerr solution is the space-time under consideration. 

These quality factors arc obtained from a number of coordinate-free characterisations of the 
Kerr solutions. A coordinate-free characterisation of the Kerr solution is a set of geometric 
conditions on a space-time which are fulfilled if and only if the space-time is locally isometric to 
the Kerr solution. A number of them have been presented in the literature, and we rely specifically 
on those which are formulated as local conditions, computable at any given point. An example 
of these can be found in [7] where the Kerr metric is characterised in terms of a set of conditions 
on the concomitants of the curvature tensor at a point, so they are fully invariant. Another local 
invariant characterisation of the Kerr solution was put forward in [2JJ . In this case one assumes 
the existence of a Killing vector in the space-time (an invariant property) and that the space-time 
is a vacuum solution. The characterization is then written as a set of local conditions involving 
the Killing vector and the curvature tensor. This result is the starting point of our work. The 
plan is to cast this invariant characterisation as a single scalar condition having the properties of 
a quality factor. Indeed, it turns out that this can be done in more than one way and we present 
a number of quality factors which are relevant in different situations. 

An important property of our approach is that the quality factors that we introduce depend 
exclusively on invariant geometric properties of the space-time under study. One does not need 
to know anything about the Kerr solution, nor one has to put the space-time in relation with the 
Kerr solution (by diffeomorphisms, for instance). Actually, the quality factors can be explicitly 
computed for any given stationary space-time. 

The question of whether or not our quality factors are good enough to measure the deviation of 
solutions from the Kerr metric will have to be evaluated, and a calibration of the meaning of any 
quality factor having values close to one (say 0.95 or more) should be carefully and exhaustively 
performed. In this sense, we support our analysis with the study of some practical examples where 
the quality factors have been explicitly computed. The results arc in agreement with intuitive 
expectations of when a space-time can be considered closed to the Kerr solution. 

The plan of this work is as follows: in section [5] we recall some generalities of vacuum space- 
times possessing a Killing vector. In sections [3] and 2] we construct a quality factor valid for 
any stationary, Ricci-flat and asymptotically flat space-time. This quality factor is based on the 
so-called Mars-Simon tensor. The more general case of stationary space-times is addressed in 
section [5] where we construct a number of quality factors based on the spacetime Simon tensor. 
We discuss quality factors valid for vacuum solutions and a quality factor which applies to non- 
vacuum solutions (although in this case the quality factor being unity is just a necessary condition 
for the space-time to be locally isometric to the Kerr solution). In section [51 we carried out 
numerical computations of the quality factors in practical cases with the aid of the system xAct 
[26l [27] . Finally we discuss some open issues of our analysis in section [Jj 
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2 Geometric preliminaries 



Let M. be a 4-dimensional smooth manifold endowed with a smooth Lorentzian metric g a b (sig- 
nature convention (— , +, +, +)). The pair (Ai, g a b) will be referred to as a Lorentzian manifold or 
spacetime. Small Latin letters a, 6, c, . . . will be used to denote abstract indices for tensors fields 
on any tensor bundle built from T(A4) and its dual T*(M). Square (rounded) brackets enclosing 
a group of indices will denote the antisymmetrization (symmctrization) operation. The inverse 
of g a b is g ab and these quantities can be used to raise and lower tensor indices in the standard 
way. The Lcvi-Civita covariant derivative compatible with g a b is the operator V a , while R abc d 
and R a b denote respectively the Riemann and Ricci tensors arising from V a . Our conventions in 
the definitions of these quantities are set by the relations, 

V a VbZ c - Vb^ a Z c = R bad c Z d , R ab = R acb c , 

where the first relation is just the Ricci identity valid for any vector field Z c . Sometimes we will 
denote tensorial quantities with index-free notation in which case boldface capital letters will be 
used (the rank of the tensorial quantity should then be clear from the context). The complex 
conjugate of a complex quantity will be denoted with an overbar. 

A vector field £ on Ai is called a Killing vector if it fulfills the differential condition 

£^g a b = 0. (1) 

where £^ denotes the Lie derivative with respect to £. The previous expression can be written 
in terms of the Killing 1-form ( a = g a bC b m the well-known equivalent form 

V Q Cb+V 6 Ca=0. (2) 

From now on we will assume that the spacetime on which we work possesses a Killing vector. 
Killing vectors have been extensively studied and their properties can be found in many places of 
the literature but in order to make our presentation as self-contained as possible we will review 
here those properties needed in our work. First of all we assume that the Killing vector £ is a 
smooth vector field in M. (indeed, the smoothness of the Killing vector follows if one just assumes 
that it is a C 2 vector field, see below). The condition © implies that Ca has the property 

V[ Q Ci>] = VaCb: 

which enables us to define an exact 2-form, called the Killing 2-form, by F a b = S7 a Cb- It follows 
immediately that 

V [a F 6c] = 0. (3) 

Using the Killing condition @ and the Ricci identity one may compute the covariant derivative 
of the Killing 2-form, yielding 

V b F ac = -R ac bd( d - (4) 

The norm of the Killing vector is £ a £ a and no restriction on it is imposed at this stage. If we 
set A = £ a £° then the Killing condition @ easily gives 

V b A = 2C a V 6 Ca = 2F ba ( a . (5) 

The identities discussed so far hold regardless of the space-time dimension and its matter content. 
In the particular case of Ricci-flat 4-dimcnsional space-time we have additional properties to be 
detailed next. First of all, in a Ricci-flat space-time, that is with vanishing Ricci tensor R ab = 0, 
the Riemann tensor coincides with the Weyl tensor C abc d which in four dimensions has the 
algebraic properties C* bcd = *C a bcd where the Weyl tensor right and left duals are defined 
respectively by 

'-'abed = "^Vcdpq'-'ab i '-'abed = "^Vabpq^ ab j 

and rj abc d is the canonical volume 4-form. Therefore, we will only use the Weyl tensor and its 
right dual from now on. Also, in four dimensions one can define the dual of the Killing 2-form F ab 
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by F* b = F pq r) abpq /2. From F* b we introduce the twist 1-form cu a = F* b ( b which, for a Ricci-flat 
spacetime, can be shown to be a closed 1-form, V[ Q Wh] = 0. The local potential arising from uj a 
is called the twist potential us. The 2-form F* b has the differential properties 

Vc^ft = C* abcd ( d => V [a F b * c] = 0. (6) 

It is possible to achieve a great simplification in many computations by introducing complex 
quantities. In this way we define the self-dual Wcyl tensor, and the self-dual Killing 2-form by 
the expressions 

Cabcd = Cabcd + i C* bcd , F ab = Fab + '^ab- v ) 

These quantities fulfill the well-known algebraic properties 



Cabpq — ^VpqedCab — i C ab p q , F ab — ^Tj abcd F — 1 F ab (8) 



and the identities 



F ac F b c — -(F ■ F)g a b , F ac Fb c — Fb c F a c , F ■ T = T a bF ah . (9) 

With the aid of these new quantities, one can write the differential conditions fulfilled by F a f, and 
F* b in a compact form 

V c .F o6 = -C abcd C, d , V [c F ab] = 0. (10) 
Also, if we introduce the Ernst 1-form 

a a = 2T ab ( b , (a a C = 0) 

one has that in a Ricci-flat space-time 

V [a< r b] = , (11) 

which means that one can find a local potential a for o~ a , called the Ernst potential, which can 
be defined in terms of the Killing norm and twist as 

(7 = l + A + 2ia;=> V a er = a a . (12) 

More generally, the Ernst 1-form is known to be closed in any space-time that satisfies C[cRa]bC b = 
0. The reader is warned about the non-standard choice of the real part of a in (|T12j) . as this 
is universally chosen to be A in the literature. However, our simple choice will simplify the 
appearance of many expressions below. Observe that a is defined up to the addition of an 
additive purely imaginary constant, so that we have the gauge freedom 

a — >a' = a + ia, ael. (13) 

This freedom will be relevant later when we deal with spacetimcs locally, without reference to 
their possible asymptotic properties. 

The quantities F a b and o~ a can be related by the formula 

- XF ab = ( [a a b] + ^rj abcd ( c a d . (14) 

See [131 123| for a more complete review of the properties of the quantities just introduced. 

3 The Mars-Simon tensors and their properties 

At those points of AA in which a ^ wc can introduce the following rank-4 tensor 

Sabcd = Cabcd + — (,F ab F c d F • Fj , (15) 
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whore 

X-\ d =^(ir, ab cd + S a c S b d -5 a d S b c ) (16) 

is the identity tensor in the space of self-dual 2-forms. Observe that S a bcd is affected by the 
freedom (TT3"| so that one has to either (i) try to choose a preferred Ernst potential, or (ii) add 
the arbitrary constant a and consider one-parameter families of potentials o~(a) and of tensors 

S a bcd{0L). 

The tensor S a bcd was first introduced in [14J 1 ! and it was called the Mars-Simon tensor, a 
terminology which we also adopt here. This tensor satisfies a number of remarkable algebraic and 
differential identities whose proof can be found in [TT] [TT] . The most relevant ones for us come 
from the fact that the Mars-Simon tensor is a Weyl field (also called a Weyl candidate), which 
means that it fulfills the same algebraic properties as the Weyl tensor 

Sabcd = —Sbacd , Sabcd = S c dab j S abc d + S a cdb + S a dbc = , S a bac =0 (17) 

and it is self-dual 

1 Sabcd = 1 *S a bcd = S a bcd- (18) 

Its covariant divergence is given in Ricci-flat spacetimes by [TTJ [TT] 

aV a S a bcd = 2C(2(T p [d S c]bap + F b p S dcap ) + g 6[ d<Sc]a P9 ^ P9 )- (19) 

We note the important property that this covariant divergence is a linear expression in the Mars- 
Simon tensor itself. Using this one can further show that S a bcd fulfills a hyperbolic equation |14| 
and that it vanishes in a globally hyperbolic space-time if it is zero on a Cauchy hypersurface 
(causal propagation) [TT] . 

The tensor S a bcd plays an important role in characterizing invariantly the Kerr solution. One of 
such characterizations requires that the space-time be asymptotically flat as originally formulated 
by Mars in [2TJ, Theorem[JJbelow. To that purpose, we briefly introduce the concept of stationary 
asymptotically flat end. 

3.1 Asymptotically flat end .Moo 

Assume that there exists a submanifold .Moo C M diffcomorphic to M. x {x £ K 3 : \x\ > R} for 
some R large enough. Supposse further, that in the local Cartesian coordinates (t, x l ), i = 1,2,3 
onAioo defined by the diffcomorphism we have C = dt (time translation at infinity) and also the 
conditions 

\9ab - Vab\ + \rd xi (g ab )\ < Ar' 3 , d t g a b = 0, (20) 

where A and B are positive constants, g ab and rj ab represent, respectively, the components of the 
metric and the Minkowski metric in the above local coordinates and r = \J (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . 
A set Moo with these properties is called a stationary asymptotically flat end. It is well-known 
(see e.g. [5]) that, using the vacuum Einstein equations (without cosmological constant), the 
components of the metric tensor have the following asymptotic behaviour in .Moo 

2M AS i x k 
gtt = -l + ^r + 0(r~ 2 ) , g txi =-e ijk — — + 0(r" 3 ) , g x * xj = ^ + O^" 1 ), (21) 

where M is the Komar mass [TS] of the Killing vector £ in .Moo, £yfc is the Levi-Civita symbol 
in dimension 3, and Sij is the Kronecker symbol (Riemannian flat metric) in dimension 3. 

3.2 Two theorems on S a b c d 

We are now ready to present the first Mars' theorem. 

lr The definition of 1141 uses a different choice of the Ernst potential. 
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Theorem 1. Let (M.,g a b) be a smooth Ricci-flat spacetime possessing a Killing vector £ and a 
stationary asymptotically flat end Moo such that £ tends to a time translation at infinity and the 
Komar mass of C at does not vanish. Fix the freedom H3\) by letting a — > at .Moo and 

construct the tensor S a bcd according to U5\) . 

If <S abed vanishes, then (A4,g a b) is locally isometric to the Kerr spacetime. 

This theorem includes a global condition, namely, the existence of a stationary asymptotically 
flat end .Moo- This is important to fix the freedom (| 13[) for the potential a to be used in the 
definition of S a b c d, as one can choose the potential such that a goes to at the asymptotically flat 
end fixing the constant a in p3[) so that Theorem [1] possesses a totally unambiguous meaning. 

Nevertheles, it is possible to formulate a purely local characterization of the Kerr metric, see 
theorem 1 in (53] . To that end, we need to control the change of S a bcd under (fT3")) , and thus we 
introduce the tensor 

Qabcd = 6 ( T ah Ted - ■ A . (22) 



The algebraic properties of T a b and T a bcd guarantee that Q a bcd is a Weyl candidate too. It is a 
very simple matter to check that under the freedom (|I3[) the Mars-Simon tensor changes as 

i icy 

S abed = Sabcd ', ~. — 7 Qabcd (23) 

0(0 + ia) 

Therefore the existence of a potential such that the Mars- Simon tensor vanishes is equivalent to 
the proportionality of S a bcd — constructed from any given a ^ — and Qabcd, with some given 
properties of the proportionality function. This will be relevant in section [5] and allowed Mars to 
prove his second, purely local, theorem. 

Theorem 2. Let (M,g a b) be a smooth non-flat Ricci-flat spacetime possessing a Killing vector 
£ and assume further that the tensor S a bcd is well defined on M.. If there exists a potential a 1 in 
such that the following conditions are fulfilled on Ai 

S'abcd = 0, T-T=-^, M6l\{0}, (24) 
then Ai is locally isometric to the Kerr space-time. 

Remark 1. The constant M here is the Komar mass of £ as introduced in subsection 13. II Note 
that the imposed condition M ^ on the Komar mass excludes the flat Minkowski solution even 
though strictly speaking it belongs to the Kerr family. In this work "Kerr space-time" will mean 
any member of the Kerr family with non-zero mass parameter. 

3.3 The space-time Simon tensor 

The vanishing of S' a bcd is of course equivalent to a particular proportionality between C a bcd and 

Qabcd- 

S'abcd =0 Cabcd = 7 Qabcd 3o G K I C a bcd = T^Qabcd (25) 

a' a + ia 

and this is why it will be convenient, and very useful, to use another tensor which measures such 
a proportionality. To that end, let us define 

S a bc = la[bCc]mrd( m J :rd + ^C™ C Cmar[cO~b\ (26) 

where here and in the rest of the paper we use the abbreviation 

lac = CaCc - ^9ac ■ (27) 

The tensor S a bc is called the space-time Simon tensor |21| and it can be defined in any 
spacetime with a Killing vector, independently of the existence of a potential a for a a - This will 
be very important later. S a bc has the algebraic properties of a Lanczos potential 

S b bc = i $a[bc} — S a bc , S a bc + Sbca + S ca b = (28) 
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hence the space-time Simon tensor is a Lanczos field (or a Lanczos candidate). From (1261) we 
deduce the additional properties 

CS abc = , ( b S abc = (29) 

so that S a bc is totally orthogonal to the Killing vector £ from which it is constructed. 

A fundamental property of this tensor is (see considerations leading to eq. (16) in [2TjV 

Lemma 1. In arbitrary spacetimes with a non-null Killing vector £ such that J- ■ J- 0, the van- 
ishing of the space-time Simon tensor is equivalent to the proportionality of the self-dual tensors 

Cabcd and Qabcd'- 

Sabc = Cabcd = LQ a bcd 

where L is a complex function. 

As seen above in (|25[) . the vanishing of S' a bcd is equivalent to the proportionality between C abc d 
and Qabcd whenever there exists a potential a for the Ernst one-form — so that the Mars-Simon 
tensor can be defined — . Thus, Lemma [1] implies the existence of a precise relation between the 
tensors S abc and S abc d- This relation must be invariant under the freedom (|13[) . because so is 
Sabc, and given that S a bcd changes according to (|23[) this relation will not involve Qabcd- This 
leads to exactly the same relation (|26[) but simply replacing C a bcd by S a bcd- 

Sabc = la[bS^mrdC m J :rd + 4C" C Smar[c°~b] (30) 

The main relevance of the space-time Simon tensor resides in the following result (lemma 4 
in m) 

Proposition 1. Let (A4,g a b) be a Ricci-flat, non-locally-flat, C 3 spacetime possessing a Killing 
vector £ which is non-null on a dense subset of M. and assume that the corresponding spacetime 
Simon tensor vanishes everywhere. Assume also that T ■ T 7^ everywhere. Then, the Ernst 
1-form is exact a a = V a <7 and furthermore: 

Cabcd = — r— Qabcd, T ■ T = b(a + a) 4 

a + a 

where a and b are complex constants. 

Combining this result with Theorem [2] we obtain another version of the local characterization 
of the Kerr space-time. 

Theorem 3. Under the same hypothesis of Proposition^ (Ai,g a b) is locally isometric to the 
Kerr space-time if and only if the constant a is purely imaginary and b is real and negative. 



4 The super energy of S ahcd 

Given a real Weyl candidate W a bcd, we define its Bel- Robinson tensor T a b c d{W} by 

Tabcd{W} = W a p d m W bpcm + W a m c p W bmdp - \ 9 abgcdW mnpl W mnpl . (31) 

o 

As described in the general framework of |32| , this tensor is the basic superenergy tensor of the 
Weyl candidate W a bcd in 4 dimensions. We prefer to call this tensor the Bel-Robinson tensor of 
Wabcd in analogy with the Bel- Robinson tensor constructed out of the Weyl tensor [3J [32] . 

In the present context it is more convenient to re-write pip in terms of the self-dual tensor 
yVabcd = Wabcd + ^W* bcd corresponding to the Weyl candidate and its dual. In this way one can 
prove that the tensor ((3T|) can be rewritten simply as (see e.g. [28] ) 

Tabcd{W} = W a P c m Wb P dm = Tabcd{W}. (32) 

As any other superenergy tensor, the Bel-Robinson tensor of a Weyl candidate has the following 
properties — see [35] for detailed proofs. 
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Theorem 4. If T abmn {W} is the Bel-Robinson tensor of the Weyl candidate W abmn then 
1- T [abmn) {W} = T abmn {W}, T a amn {W} = 0. 

2. Generalized dominant property: if u\, u*2, ua are causal future- directed vectors then 
Tabmn{W}u1u2U™u' l l > 0. The inequality is strict if u\, H2, U3, ua are timelike. 

3- T abmn {W} = <^=> Wabmn = <^> W abmn = <^> 3 a set of timelike vectors ui, u 2 , 
M3 and H4 such that T abmn {W}uiU2U™u'l = <S=> there is a timelike vector u such that 
T abmn {W}u a u b u m u n = 0. 

The last point of this theorem enables us to write the vanishing of the self-dual Weyl candidate 
yVabcd as a single scalar condition. For an arbitrary timelike vector u we define the scalar 

U S (W) = T abmn {W}u a u b u m u n (33) 

which, if Wabcd 7^ 0, is a positive quantity according to points[5]and[3]of Theorem^ Therefore one 
has a natural scalar quantity which enables us to measure, for each timelike vector, the proximity 
to the geometric conditions determined by the tensor condition W a bcd = 0. The scalar quantity 
Ua{W) is called the superenergy density of T abcc i{W} with respect to u }32pl and it depends in 
general on the timelike vector u. 

However, in those cases where there is a timelike vector defined invariantly one can select a 
superenergy density that acquires an invariant meaning. We can apply this idea to the particular 
case in which W a bcd = S a bcd- At those points of M. in which the Killing vector £ is timelike 
the scalar Ug(S) is positive and it will vanish if and only if S a b c d = 0. Given that £ is defined 
invariantly we can take the quantity U^(S) as a local invariant measure of the deviation of the 
spacetime (Ai,g a b) to the geometric conditions entailed by S a b c d = 0. 

Theorem 5. 

UJS) = £ ab £ ab + £j|p (3(<r • <x) 2 - \cr ■ ct\ 2 + 4Rc(a£ ab a a a b )) , (34) 

where 

£ab = C ap b q ( P ( q , <T • <T = (J a (j a , <T • & = 0~ a (j a . (35) 



Proof. Starting from (|15[) and using the definition of the Ernst 1-form, a straightforward compu- 
tation yields 

U ac = S abcd (\ d = £ ac + ^7ac + (36) 

Za la 

Next, we write the definition of U^(S) adapting the definition of T abc d{S} to eq. (|32| . The result 
is 

U c (S)^U ac U ac . (37) 

Now we need to replace in this expression the value of U a b found above and work out the resulting 
expression. After a computation, we find 

U ac U ac = £ab£ ab +^£abO- a a b +^£ ab a a a b + -^ ( iia-af-XiF-Fcr-a+F-Fa-^+X^F-FA . 

la la % tJ \ \ / 

(38) 

The final expression (|34p follows after using in this expression the identity 

A T- T = cr cr, (39) 
which is a consequence of (|14|) . □ 
Remark 2. Note that the result is valid for a Killing vector £ of any causal character. 

2 The notation differs from that of 1321 . 
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One may also use the previous definitions and results to introduce the quantity 

U^Q) = T abcd {Q}CC b CC d (40) 

for the Wcyl candidate ((22|) . According to Theorem [H the scalar Ug(Q) will be non-negative 
whenever £ is causal. Indeed a straightforward computation yields 

U^Q) = ^-^^^\\ (41) 
Hence (j3"4"]l can be rendered in the form 

UW-UwM + '^'ff* . (42, 

<• C |of |cr| 2 

4.1 Several alternative expressions for the scalar Ug(S) 

The superenergy construction can be carried out for arbitrary tensors of any rank, and not only 
for the Weyl candidates [32]. For a 2-form F ab the corresponding superenergy tensor T ab {F} is 
simply the standard energy-momentum tensor of a Maxwell field and can be expressed in any one 
of the following equivalent forms 

T ab {F} ee t ab = F ac F b c - ±g ab F cd F cd = 1 (F ac F b c + F* C F^ C ) = ^F ac F b c . (43) 

Using these equations and (0) we can easily derive the identities 

tab = ha , t a a — , t ac t b c — -g a b t cd t cd , t cd t cd = — ( F ■ F) ( T ■ F) , 

4 lb 

and by virtue of the "Maxwell" equations V a F ab — 0, t ab is divergence free 

V a U = . (44) 

As any other superenergy tensor, it also possesses the dominant property, that is to say, t ab u1u2 > 
for arbitrary future-pointing vectors u\ and ui , and for timelike u one has 

t ab u a u b > 0, t ab u a u b = tab = <f=> F ab = . 

On using Ea.(|14|l the tensor t ab can be written in our case as 

tab = [2(<r • <r)CaC& + A (cr o 0- 6 + a a a b ) - A(er ■ a-)g ab + 2C( a «i>)] (45) 

where we have introduced the real one-form 

v b EE ii] bcde a c ( d a e 

which obviously satisfies 

v b ( b = v b a b = v b a b = 0, v b v b = A [(er ■ <r)(<x ■ ct) - (cr ■ ct) 2 ] . 

With these definitions at hand, we can find an interesting expression for £ ac . To that end, 
simply use the first in (fTU)) to compute 

£ ac = -C b V c F ab = -\v c „ a + W = -\v c a a + \F ab (F b + F b ) 

so that introducing here the identity © together with (|39p and the definition of t ab we arrive at 
the desired expression 

£ ac = t ac + -^r(cr ■ cr)g ac - -V c cr a . (46) 
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Observe that this provides a very simple formula for the electric and magnetic parts of the 
Weyl tensor with respect to the preferred timclikc direction £ in Ricci-flat stationary spacctimes 
in terms of the Ernst one-form. 

Using (f2"T)l and we can also derive 

1 , x 1„ (c ' 3(7 a <7 c .,„. 

% c = i ac + — (a- ■ a)g ac - - V c cr a + ^Aac + . (47) 

From the previous formula (|46[) it is very easy to compute the Bel-Robinson superenergy with 
respect to the direction £: 



UAC) = - [V a * b V a a b - 2t ab V a (a b + a b )] = - [V a V b aV a V b a - 4UV a V b A] 



which on using (|45[) can also be written in terms of the Ernst one-form exclusively as 

U Z {C) = ^W- I ^(,r. ( T)(a.a) + ^ [{ct ■*) + (*■ a)]-±a a a b V a (<T b +<T b ). (48) 

This is an expression for the Bel-Robinson super-energy relative to the direction £ involving only 
the Ernst one-form and its derivatives. 

We can also use together with (|4U)) to compute 

£ ab a a a b = [{a ■ a) + {a ■ a)} - \a a a h ^ b a a . 



The previous two expressions together with (|42Mip lead after some manipulations to 
Ut(S) = ±V^W + %^ - (6A* - 6A + W\ 2 )^ + 
2Re 



; (6A + <t)ct - ^^V b( x a - — V 6( r a 



32A 2 |a| 2 v 8A " " 4a 



(49) 



which again is an expression depending exclusively on the Ernst potential and its derivatives. 
We can summarize the previous results in the following theorem 

Theorem 6. Let [M., g ab ) be a smooth spacetime with a Killing vector £ that contains a stationary 
asymptotically flat J^-end M.^ such that C, tends to a time translation at infinity and the Komar 
mass of C, at Moo does not vanish. Fix the freedom U3\) such that a — > at Moo. 

Then, the scalar J^ff| ) constructed exclusively from a and its first and second derivatives is 
non-negative on the entire region where C * s timelike, and it vanishes if and only if (M,g ab ) is 
locally isometric to the Kerr spacetime. 



5 Quality factors that measure the Kerr-ness of stationary 
spacetimes 

In this section we present several quality factors to measure how close a given spacetime is to the 
exterior part of the Kerr solution. By a quality factor we mean an adimcnsional scalar function 
q with q G [0, 1] and such that its value provides a measure of the "Kerr-ness" of the given 
space-time at any point. 

We present three different possibilities. All of them can be applied to generic stationary 
spacetimes, the first two to stationary Ricci-flat ones, while the first one can only be used when 
the space-time is Ricci-flat and has an asymptotically flat end. 
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5.1 A quality factor for stationary Ricci-fiat asymptotically flat space- 
times 

Theorem [6] provides a neat characterization for the Kerr solution in terms of a single scalar and 
real quantity, which happens to be strictly positive if the metric is not Kerr. One can thus think 
that the scalar (|49p may provide a measure of how much a region of the spacetime differs from 
the stationary portion of Kerr spacetime: the smaller the scalar is, the better resemblance 
to the Kerr solution. Unfortunately, the previous sentence does not have an absolute meaning, 
because the scalar (|49| has physical units — of (length) -4 — and thus we can always choose the 
physical dimensions such that it becomes large or small. 

What is needed here is another similar quantity that normalizes (|49[) and eliminates the 
physical units. There may be several choices here, but our own, which we believe is a natural 
one, is based again on the superenergy construction and its positivity properties. As a matter 
of fact, one can define a (positive-definite and sesquilinear) inner product ( , )j on the complex 
vector space of self-dual Weyl candidates as follows: 

(Wi, w 2 ) e = Wi a V" W 2 bpdm cc b CC d 

Scsquilincarity follows obviously from the definition, while its positive definitcness is a direct con- 
sequence of Theorem 01 Observe that the norm here is simply the superenergy scalar introduced 
in (23]) 

\\W\\l=(W,W) i = U i (yV). (50) 

Definition 1 (Quality factor q). The quality factor measuring the Kerr-ness of any stationary 
region in a non-flat asymptotically flat and Ricci-fiat spacetime (M,g a b) is defined by 



where 

Q 2 



q=l-Q Z 

UAS) 



uac) + \o-\-^Ju~ c {q) 



The main properties of this quality factor are collected in the following 
Theorem 7. The quality factor q has the following properties: 

1. ?e [o,i]. 

2. q — 1 on an open neighborhood U C M if and only if {U, g a b) is locally isometric to the 
Kerr solution. 

3. At any point x £ M, q~ if and only ifC a bcd\x = k{Qabcd/o')\x for a positive real constant 
k € R+. 



Remark 3. Notice that the condition for q = is that the space-time is of Petrov type D — as it 
is the case for the Kerr solution — , however the proportionality constant k has the opposite sign 
to the case of Kerr. 

Proof. From ([50]) we can rewrite Q 2 as 

II S ||| 

Q 2 = - c - — - 



(lie || C - + 11 1 



"c 



and using here that (indices suppressed) S = C + Q/a the triangle inequality 
\\C + Q/o-\\z<\\C\\£+\\ Q/a ||j leads to 

< Q 2 < 1 
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which proves 1. Now, Theorem |B] informs us that the space-time is locally isometric to Kerr if 
and only if Ug(S) =\\ S ||S= 0, that is, if Q 2 — 0, which proves 2. To prove 3, recall that the 

triangle inequality is an equality if and only if Q/<j\ x and C\ x are positively linearly dependent, 
i.e. if C\ x = k{Q/o)\ x with k > 0. The case with k — is excluded as the space-time would be 
flat. □ 

Remark 4. It is very important to remark that q can be computed in any given non-flat Ricci- 
flat and asymptotically flat stationary space-time without any reference to the Kerr solution. One 
does not need to map the space-time into the Kerr solution, nor is it necessary to compare or 
choose coordinates relative to Kerr. All that one needs is an explicit expression of the metric of 
the space-time to be compared to the Kerr solution. Notice, furthermore, that using (|49 p .(|48 p 
and (|4ip one can immediately write an expression for q in terms of a and its derivatives. 

At any point x of the stationary region in (A4,g a b) one can say that x is "(g x 100)% Kerr" 
(so it is 100% Kerr if q = 1, and 0% Kerr if q = 0). The precise value of q which gives a good 
approximation to Kerr depends on the context, and it will have to be calibrated by studying 
several examples. Some of these are given in Section [BJ 



5.2 Quality factors for stationary Ricci-fiat spacetimes 



There may arise situations where one does not know if the entire space-time is asymptotically flat, 
or also when one knows that there are no asymptotically flat ends but still the proximity with a 
Kerr space-time is suspected, for example close to the ergo-surface of the Killing vector or to a 
possible null hypersurfacc which one suspects may resemble the horizon of Kerr space-time. In 
these situations we cannot apply the previous quality factor, which requires fixing the constant a 
in ([13")) for the potential a. However, we can still make use of the previous results by complicating 
things a little bit, as follows. 

For Ricci-fiat spacetimes we can make use of Theorems[2]and[3]which do not require asymptotic 
flatness. The first thing we must require to use those theorems, according to (|2"5)l or Lemma [TJ 
is the proportionality between C a bcd and Q a bcd- Note, however, that this can be easily measured 
by means of 



(C, Qh 



eii| 



Theorem 8. The quality factor q\ has the following properties 

1. qi e [0,1]. 

2. q\ = 1 at a point x € M. if and only if 31 G C : C a bcd\x = IQabcdl 

3. qi = at a point x E M if and only if £ a ba a a b \ x = 0. 



Proof. The Cauchy-Schwarz inequality immediately gives 

< qi < 1, qi\ x = 1 <^=> Cabcd\x = lQabcd\x, (51) 

for some I € C. The other extreme value, qi = 0, can only occur when 

(C, Q) c - - \E a bO a o h 

vanishes, proving point 3. □ 
Point 2 in this theorem informs us that q\ measures how much C a bcd and Q a bcd are proportional 
to each other. A possible alternative to measure this proportionality is the use of the space-time 
Simon tensor S a b c - We first note the following simple result. 

Proposition 2. Assume that the Killing vector C, is time-like. Then 

SabcS ahC > , S a bcS abC = S a bc = 0. (52) 
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Proof. Since C is time-like, from we know that iS a f, c is totally orthogonal to £, and thus 
spatial with respect to the Killing vector, from where the result follows easily. □ 
This result together with Lemma [T] enables us to use the quantity S a b c S a in a definition of 
the quality factor to measure the proportionality of C a bcd and Q a bcd- To that end, for any two 
Lanczos candidates S\ bc and S 2 bc (not necessarily traceless), define 

(S\S*)^Sl bc S 2a " c . (53) 

This is again a sesquiliear product on the space of Lanczos candidates (not necessarily traceless), 
and it is positive definite when acting on candidates that are totally orthogonal to a timelike 
Killing vector £, as follows from Proposition [2] Actually, 

ll-sil 2 ^ (S,S) 

happens to be essentially the superenergy of the (real and imaginary parts) of S a bc with respect 
to Ci according to the general construction in [3"2] . 

Now, observe that S a bc can be written, from its definition ((26]) . directly as the sum S\ bc + S 2 bc 
where 

s lbc = --jla[b^c\ ^ ' ^ ' S * bc = _4 ^»[i> f7 c]> 

which are Lanczos candidates on their own, and they are totally orthogonal to £■ Therefore we 
can define 

92Sl II*" 2 ,. (54, 

(II S> || + || S' II) 1 

Theorem 9. The quality factor qi has the following properties 

i- 92 e [0,1]. 

2. q2 = 1 at a point x e M if and only if3l€C: C a bcd\x = lQabcd\x- 

3. qi — at a point x £ M. if and only if V c (J- ■ J-) \ x = (thus J- ■ J- = const, if qi =0 on 
an open neighborhood) which entails S^ bc \ x = 0. 

Proof. The first property is a direct consequence of the triangle inequality 

|| S || = || S 1 + S 2 \\<\\ S 1 || + || S 2 ||, while q 2 \ x = 1 happens if and only if S a b c \x = as is clear 
from (|54p. so that Lemma [1] entails then C a bcd\x = lQabcd\x for some 1 6 C. To prove the third 
point, recall that q%\ x = if and only if <S* 6c | x and S 2 bc \ x are positively proportional to each 
other, that is 

Slbc\x = kS 2 bc \x, fceM, k>0 

which reads 

-jla[b^c] {F ■ J 7 ) \x = 4:k£ a[b cr c ] \ x . 
Contracting here b and a we derive 

AV C {F-F)\ x =Ak£ ac a a \ x 
but using here (|46|) and (|45j) together with (|39|) we also have 

£ac<J a = + * c ) ~ \o- h ^cO h = ^jf^ c \ - jV c (<7 • a) = -^V c (T ■ F) (55) 

so that the previous condition becomes 

(l + fc)AV c (T- T)\ x = 

which proves the result. □ 
Actually, ea. (|55[) allows us to derive an alternative useful expression for iS a f, c . 
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Lemma 2. The space-time Simon tensor can be written in the form 

Sabc = j(4:£a[blc]d + 2£d[clb]a)(? d - (56) 

Proof. By noting that r ) c d,o~ d = — A<r c the first summand provides directly S 2 bc while the second 
summand is S\ bc as follows from (|55|) . □ 
If we combine eq. (fBTj) with the second point of Theorem [§] we deduce that 

q x = 1 g 2 = 1 . (57) 

One can find an alternative, useful, expression for qi as follows. By explicit computation one can 
easily obtain the following results 

IIS 1 II 2 = jV c (J.J)V c (J-J) (58) 
<S\S 2 ) = -|| SMI 2 (59) 
||S 2 || 2 = 8 [{a ■ &)£ ab £ ab - £ b c £ ac o a a b ] (60) 

where we have used (|55p and ([551) to get (|50[) . From these one immediately finds 

(S\S) = Q, || S || 2 =|| S 2 f - || S 1 || 2 (61) 

which also implies 

||S 2 H^HS 1 || 2 . 

Observe that the second in (|6ip allows us to rewrite qi as given in (|54[) in a simpler form 

_ 1 ||S 2 || 2 -|| SMI 2 _ || S 2 || - || S 1 || 2 us 1 1| 

92 (|| S 1 || + || S^ ||) 2 || S 1 || + || S^ || || S 1 || + || S^ || 

from where the properties in Theorem [5] become more transparent. 

Remark 5. We note that the quality factors q\ and q2 do not involve the Ernst potential and 
hence they can be computed for any space-time possessing a Killing vector, be it Ricci flat or not. 

Assume then that any of q\ or q2 were unity. In that case, from Proposition [T] we know that, 
for Ricci-flat spacetimes, there would exist a complex constant a such that L above would be 
L = — 1/(<7 + a) for a potential a with Rc(er) = 1 + A (redefining a if necessary). Then we would 
have 

= -^H fill? ^ 
so that the Kerr condition a = —a (see Theorem [3]) would require that the quantity 



/3 = Re 



Qll 2 ; 



(C, Qh 



equals — (1 + A). Notice once again that to compute f3 we only need to know the space-time and 
the Killing vector, and therefore we can control this property by means of the partial quality 
f&ctor 

>6s2 l^ + lllA|-^ + l)A = + + + 

(l/3 + l| + |A|) 2 (l/^ + ll-A) 2 

where in the second equality we have used that A < 0. It is easy to prove that 
Proposition 3. The quantity k fulfills the following properties 

• < k < I, 

• k = if and only if (3 + 1 < , 
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• K = 1 if and only if j3 + 1 = —A (this is the Kerr condition). 



Proof. An alternative expression for k is 

([3+1 + A) 2 



K= 1 - 



(l£ + l| + |A|) 2 

from where the first and third properties are obvious. The second follows directly from □ 

Provided that gi = 1 (or cquivalcntly qi — 1) we must finally control the second condition on 
Theorem [31 that is to say that the constant b in Proposition [T] be real and negative. From (p2")) 
and (|2H]) this can be rewritten as 



C||§ 



being real and positive (as this would have to be z = bX). Then set 

"-'-Kh-'Xr- 1 

Proposition 4. The following statements hold true 

• < x < 1, 

• x = if and only if z = —\z\ = z (i.e. z is real and negative/ 

• x = 1/2 if and only if z — —z (i.e. z is purely imaginary). 

• x = 1 if and only if z — \z\ — z (i.e. z is real and positive/ This is the Kerr condition. 

Proof. If we set z/\z\ = e 1 ^ then elementary manipulations show that the scalar x takes the 
following form in terms of <f> 

* = cos2 ($) ' ( 64 ) 
which makes it clear that < x < 1. Now. one has 

• x — 4> = (2n — 1)tt, n G Z z is real and negative. 

• x = 1/2 <^=> <fi = tt/2 + nn, n6Z z is purely imaginary. 

• x = 1 = 2nn, n G Z z is real and positive. 

□ 

Thus, we can define two new quality factors to measure the Kerr-ness of general stationary 
Ricci-flat spacetimcs. 

Definition 2 (Quality factors q and q). Two quality factors measuring the Kerr-ness of any 
stationary region in a Ricci-flat spacetime (AA,g a b) are defined by 

q = q 2 KX, q = qiKX, 

The main properties of these quality factors are 
Theorem 10. The quality factors q and q have the following properties: 

1. q, qe [0,1]. 

2. On an open neighborhood U C M. we have that q = 1 q = 1 (U,g a b) is locally 
isometric to the Kerr solution. 

3. At any point x G Ai, q\ x — if and only if either (/? + l)^ < 0, or z\ x is real and negative, 
or V c (F -T)\ x = 0. 
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4- At any point x € M., q — if and only if either (/3 + l)| x < ; or z\ x is real and negative, 
or£ ab a a a b \ x =0. 

Proof. The first item follows from the fact that each quality factor is the product of quantities 
in the interval [0,1]. To prove the second item, we note that if q = 1, then clearly q2 = 1, k = 1 
and n = 1, as q2, K, h € [0,1]. This together with ([57)) implies that 

According to Proposition HJ Proposition [3J Theorem [S] and Theorem [HI the fact that <?2 = 1 , 
q\ = lj k = 1, x = 1 in an open neighborhood J7 C A4 is equivalent to the conditions of Theorem 
[3] being met and hence (U, g a b) is locally isometric to the Kerr solution. Similarly, points 3 and 
4 of the theorem follow directly from the mentioned propositions and theorems. □ 

Remark 6. One could also have defined q' = (q2 + K + x)/2> — and the same for q — which also 
has properties 1 and 2 of the theorem. However, observe that the chosen q is more demanding, in 
the sense that larger values of q are required to have a good resemblance to the Kerr space-time. 
For example, if tfe = 1> k = 1/2 and x = 1/2 we would have q' = 2/3 while q = 1/4. 

5.3 How to apply the quality factors to general stationary spacetimes 

The results described so far assume that (AA,g a b) is a vacuum solution of the Einstein field 
equations (without cosmological constant A) and hence they cannot be used to establish the 
proximity of a non-Ricci-flat solution to the Kerr spacetime. The problems are of two types. 
First of all, the Mars-Simon tensor can only be defined on spacetimes with a potential a for 
the Ernst one-form, and this only happens when C[bRa]cC c = 0. A way to circumvent this is to 
replace the Mars-Simon tensor with another quantity which can be defined in any spacetime, 
independently of the existence of the Ernst potential. Obviously, the space-time Simon tensor is 
an useful alternative, as it can be defined in any spacetime with a Killing vector, be it Ricci-flat 
or not, and independently of the existence of a potential a. And this is where the second difficulty 
arises, because the vanishing of the space-time Simon tensor implies the proportionality of C a bcd 
and Qabcd, but the explicit form of the proportionality factor given in Proposition [TJ as well as 
the explicit expression appearing there for T ■ J 7 , depend crucially on the assumption that the 
Ricci tensor vanishes. 

Nevertheless, it is not only desirable, but actually important, to be able to measure the prox- 
imity to Kerr for stationary spacetimes with non-vanishing Ricci tensor. And this is important 
at least in three respects: (i) to consider approximate solutions to the Einstein vacuum field 
equations and check their resemblance to the Kerr space-time, (ii) for cases with a cosmological 
constant A, and (iii) for more general cases such as solutions for a perfect-fluid energy-momentum 
tensor, or matched solutions containing a vacuum region and an interior part, and so on. In all 
cases we are going to argue that we can actually use some of the previous results conveniently 
adapted to the situation under consideration, and with some precautions. 

We split this analysis in two, the case with an Ernst potential, and the case without it. 

5.3.1 The case with an Ernst potential: R a bC = MCo 

This case includes, obviously, the vacuum case with a cosmological constant R a b = A<7 a £>, as well 
as for instance any space-time with a rigidly rotating fluid, which are characterized by having an 
energy-momentum tensor (ergo a Ricci tensor) with a timelike eigendirection proportional to a 
timelike Killing vector (. 

Now the relations of Proposition [1] no longer hold in general, as they are derived using the 
Ricci-flat condition. Nevertheless, the factors K and x can still be defined while qi and q-i still 
measure the proportionality between C a bcd and Q a bcd- One can expect intuitively that, provided 
the space-time is actually close to the Kerr solution somewhere, the deviation from Ricci-flatness 
will be small there, and the relations of Proposition [TJ will hold approximately, that is, with some 
small corrections proportional to the non-zero Ricci components. Thus, it seems clear that the 
quality factors q and q will still measure the Kerr quality of the space-time under analysis. It 
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may also happen that the space-time is asymptotically flat, in which case we can also use q. In 
all cases it should be kept in mind, however, that values of q, q and q close to 1 will be relevant 
only in combination with the smallness of these Ricci terms — compared with some characteristic 
length squared — . 

An obvious example to test the previous arguments is the Kerr-de Sitter space-time. Another 
important example is the rigidly rotating perfect-fluid solution due to Wahlquist, which happens 
to contain the Kerr spacetime as a particular case. These spacetimes are discussed in the last 
section where the quality factors are explicitly computed and plotted. 

5.3.2 The case without Ernst potential: C[bRa]c( c 

This is relevant for example for approximate solutions of the Ricci-flat equations, where despite 
the fact that the Ricci components are assumed to be small one has no extra control on them. It 
is also relevant for more general situations. 

In this case there arise some doubts about the factors n and x. One possibility to provide them 
with a reasonable meaning is to also require that the deviations from having an Ernst potential 
is very small. To than end, recall that in a general space-time the Ernst one-form satisfies [53] 

V[ Q CT fc ] = iVabcdCR d eC 

which in particular implies that V[ a <7;,] is totally orthogonal to ( b . Thus, the quantity V[ a <Tb] V^g-b] 
is non- negative and vanishes if and only if so does the exterior differential of a a . We can thus use 
the following scalar 

as an appropriate quality factor measuring the vanishing of V[ a crf,] V^a 6 ' , and thereby the local 
existence of an Ernst potential. Observe that q a = 1 if and only if there exists a such that 
o"a = V a cr locally. Then, one can use cither of the following quality factors for general stationary 
spacetimes 

Ma, qqa 

keeping in mind that only values of q a near one provide a reliable result. 

To test the behavior of this refined quality factors we have used the standard Lense-Thirring 
approximate solution for a rotating stationary isolated body and the Kerr-Newman solution of 
the Einstein-Maxwell equations, sec Section [51 

6 Examples 

In this section we present explicit numerical studies of the quality factors introduced above. 
6.1 The Curzon-Chazy solution 

The Curzon-Chazy solution is a vacuum solution to the Einstein equations which is asymptotically 
flat. It is the simplest member of the Weyl class (see [33J for further details). The explicit form 
of the metric is 

o 2M n 2M / M 2 sin 2 8 n n 0i o o o\ 

ds 2 = ~e~ — dT 2 + e~ (e ^—(dR 2 + R 2 d6 2 ) + R 2 sin 2 9dtf) , M e K (65) 

The metric functions are everywhere regular except at R = 0. The structure of this singularity 
is rather complicated, being a directional singularity, and it has been studied by a number of 
authors (see [15] and references therein for a good account of this). In the present study we 
restrict our attention to the region — oo <T<oo,0<i?, 0<6><7r and < <fr < 2ir. In 
this region there exists an assymptotically flat end A^oo located at {R > a}, with a <E K big 
enough and d/dT as the asymptotic time translation. Therefore we can use the quality factor q 
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Figure 1: Graphical representation of the quality factor q for the Curzon-Chazy solution when 
M = 12. In this picture we have defined the variables y = Rcosd, x = Rs\x\9 so the vertical line 
x = is the axis of the axial Killing vector d/d(j>. The curves are contour lines with constant 
q and the contour interval is 1/100. We can indeed obtain a contour surface by rotating each 
contour curve around the axial Killing vector axis. The white central region is a set of points 
where q < 1/100; actually, q = on the axis for R < M(= 12). On the other hand, the outer 
white region is a set of points where q > 99/100; in fact, q — > 1 when R —> oo. The plot indicates 
the presence of points whithin the axial Killing vector axis where q is not diffcrcntiable. These 
are probably due to the use of square roots in the denominator of q, nevertheless q is perfectly 
continuous everywhere. 



introduced in Definition [T] to compute the proximity of this solution to the Kerr space-time. In 
fact the explicit form of q for this solution is 

M 2 sm 2 8(e^f - l) ' (A/ 2 + 3i? 2 ) - 3A/ 3 i?sin 2 (e^ - l) + 3R 2 (e 2 f (M - R) + M + 

eTr - l) ^M 2 sin 2 9 (M 2 - 3MR + 3i? 2 ) + 3i? 2 (A/ - R) 2 + 2V3MR 



where £ = d/dT . In figurc[T]we present a graph of this quantity for a given value of the parameter 
M. Observe that q reduces at the axis of symmetry (6 = 0) to 



q\e=o = 1 



(M - R) + M + i?) 2 

(on,/ \ ^ 

\M -R\- \M -R\ + 2M 



so that q\e=o = at the axis for R < M, but q\e=o ^ at the axis for R > M. In addition, q — > 1 
when R — > oo . 
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6.2 The 5 = 2 Tomimatsu-Sato solution 



The Tomimatsu-Sato space-time [34| can be described by a metric whose non-vanishing compo- 
nents are given by 

A(x,y) 2MqC{x,y){-l + y 2 ) M 2 B(x,y) 



9tt — ~Tr, r 7 9tc/> — 



B(x,y) ' y <P B(x,y) ' * p 2 5 2 {-\ + x 2 )(x 2 - y 2 f ' 



M 2 B(x,y) 

9yy 



p 2 S 2 (-l + y 2 )(-x 2 + y 2 ) 3 ' 

M 2 (-l + y 2 )(p 2 B 2 (x, y)(-l + x 2 ) + 4q 2 S 2 C 2 (x, y)(-l + y 2 )) 



(66) 



A(x,y)B(x,y)5 2 
M , p, q, S being constants, and 

A(x, y) = (p 2 (x 2 - I) 2 + q 2 (1 - y 2 ) 2 ) 2 Ap 2 q 2 (x 2 l) (l - y 2 ) (x 2 - y 2 f , (67) 

B(x, y) = {p 2 x i + 2px 3 - 2px + q 2 y 4 - l) 2 + Aq 2 y 2 (px 3 - pxy 2 - y 2 + if , (68) 
C(x, y) = p 3 x (1 - a; 2 ) (2 (x 4 - l) + (x 2 + 3) (l - y 2 )) + 

p 2 (x 2 - 1) ((x 2 - 1) (1 - y 2 ) 4x 2 (x 2 - y 2 )) + q 2 (l - y 2 ) 3 (px + 1). (69) 

This is a vacuum solution of the Einstein equations if S = 2 and p 2 + q 2 = 1. The coordinate 
ranges are — oo < t < oo, —1 < y < 1, — oo < x < —1, 1 < x < oo. In the figure[2]a numerical 
study of the quality factor q for £ = <9/9i is shown. 



6.3 The Wahlquist perfect-fluid space-time 

A rigidly rotating perfect fluid solution to Einstein field equations was found by Wahlquist in 
[35] . see also [36l [31], which is known to have a vanishing Simon tensor [19l [31] and is the 
unique stationary and axially symmetric such solution of Petrov type D with a particular linear 
equation of state [3D] . The Wahlquist spacetime actually contains the Kerr-de Sitter space-time 
as a particular limit case [21] . There are several equivalent but different forms for the Wahlquist 
solution in the literature, see [55J[5DJ for a discussion of this, and the best choice for our purposes 
is the following in local coordinates {t, <p, y, z} 

a V-U_ ( ViV + v 2 U j V , UV(vi+V2) J a , , , ( dy 2 dz 2 \ 

ds = -^TVA v-u d(p ) + v-u dtp + ^ + V2) {— + —) 

where vi(z) = (cosh(2/3z) - l)/(2^ 2 ), v 2 {y) = (1 - cos(2/3y))/(2^ 2 ), and 

tt( \ n f , , smh(2/3z) ( 

u W = Qo-\vb + p)v 1 + —y—\a 1 + z-p) 

Trf \ n , ( , Mo\ . sin(2^y) / fi \ 

with Qq, ai, 0,2, vq, fJ-o and (3 arbitrary constants. The Kerr-de Sitter-NUT line element is obtained 
by setting j3 = 0, and a\ is essentially the NUT parameter. Thus, good values of the quality 
factors are only expected for small /3 and a\. For this solution, the space-time Simon tensor 
vanishes identically and therefore the quantity q± (and hence q 2 ) is identically 1. Also the Ernst 
1-form is closed so there is a local Ernst potential (or equivalcntly q a = 1). This means that 
qq a = qq a = kx. We provide a numerical study of nx computed for the Killing vector £ = d/dt 
in figure [3] 



6.3.1 The Kerr-de Sitter spacetime 

As said above the Kerr-de Sitter spacetime can be regarded as a particular case of the Wahlquist 
solution just discussed. The Kerr-de Sitter limit is obtained if we set the restrictions [53] 

Oo = A 2 , vq = 1- A 2 A/3 , at = , a 2 = -2M , /i = -A , p = , 
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Figure 2: Representation of the Tomimatsu-Sato global structure for p = 1/5. The grey region 
corresponds to those points in which the axial Killing vector d/d<fi is space-like, so no causality 
violations are expected there. The variable X is defined by the condition X = —1/x. The lines 
X = ±1 correspond to conical singularities and the lines y = ±1 are rotation axes for the Killing 
vector d/d(f>. The vertical line X = is space- like infinity for two diferent asymptotically flat 
regions (these are each characterised by —1 < X < and < X < 1). The two black thick 
closed curves are the ergosurfaces and the two dots are metric (naked) singularities. We have 
superimposed some contour curves for the quality factor q indicating the corresponding value of 
q on each curve. In this particular case this scalar turns out to be positive everywhere (and it 
approaches q^ = 1 at X = 0). 



where A, M and A correspond to the angular momentum per unit mass, the mass and the 
cosmological constant of the Kcrr-dc Sitter spacctimc, respectively In figure 2] we show a plot of 
the quality factor q = q = kx computed for the Killing vector £ = d/dt. 

6.4 The Kerr-Newman spacetime 

Consider the Kerr-Newman spacetime in the standard Boycr-Linquidst coordinates 

_ 2MR - A 2 cos 2 0-Q 2 -R 2 A(2MR - Q 2 ) sin 2 9 

gtt ~ R 2 + A 2 cos 2 9 ' ^ ~ R 2 + A 2 cos 2 6 ' 

, , sm 2 9((A 2 +R 2 ) 2 ~A 2 (A 2 + Q 2 + R 2 -2MR)sm 2 9) 
9ee = R 2 + A 2 cos 2 9, g„= i WTA^9 ' 

R 2 + A 2 cos 2 9 

9RR ~ A 2 + Q 2 +R 2 -2MR (70) 

which is a solution of the Einstein-Maxwell equations. For the Kerr-Newman spacetime the Ernst 
1-form is not closed and hence one cannot define the Ernst potential as used here. However, as is 
well known }33) . one can define a complex potential — in general Einstein-Maxwell spacetimes — 
for a linear combination of the twist one-form u> c with a one-form associated to the complex 
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Figure 3: In this figure we show a representation of the quality factor q = q = kx for the Wahlquist 
solution with parameters Qo = 4, a\ = 3.015, 02 = 1.585, vq = 0.64, fj,Q — 1.265, /3 = 1.88 and 
C = d/dt. Since we are only interested in the region which is close to Kerr, only contour curves 
at which q > 90/100 are represented (in this case the contour interval is 5/1000). The picture on 
the left is a representation of the factor kx in the plane z — y. To facilitate the understanding 
of this graph, the coordinate z has been compactificd by means of the transformation z — > tan z 
and the grey region is the region where £ is causal. The black vertical lines correspond to points 
where the determinant spanned by the Killing vector vanishes (note that the leftmost of these 
lines serves as the boundary of the grey region). The pattern shown in this picture repeates itself 
quasi-periodically along the y-axis. The black dots within the grey region represent points where 
q > 90/100. The graph on the right is an inset of a region containing the black dots, made to 
better appreciate the structure of such region. These numerical representations indicate that q 
reaches the unity value only at isolated points within the studied domain. 



self-dual Maxwell 2-form and the Killing vector. Using this generalized potential one can define 
generalizations of the Mars-Simon tensor and find characterizations of the Kerr-Newman space- 
time analogous to those used in this paper for the Kerr case: see [37] for the details. Given 
that the Kerr-Newman space-time is asymptotically flat one can thereby define a quality factor 
completely analogous to q. This is, however, outside the scope of this paper. 

On the other hand, we can pursue our program and compute the quality factors qq a and qq a 
for the Kerr-Newman solution. It turns out that q\ — q2 — 1 as one can explicitly check, and 
therefore, the quality factor which we need to study is qq a = qq a = q a Kx. In this case one can 
compute the factors K and x in a closed form getting 

O 4 

k=1 2, (71) 

(A 2 cos(20) + A 2 + 2R{R - 2M) + 3Q 2 ) 2 



(A 2 M cos 2 6+R(MR-Q 2 )) 2 , , g g 



_ _J (A 2 cos 2 8+R 2 )(A 2 M 2 cos 2 6+{Q 2 -MH) 2 ) 

A 2 Q 4 cos 2 8 

(A 2 cos 2 e+R 2 )(A 2 M 2 cos 2 6 + (Q 2 -MR) 2 



Si § n (9(i^))<0 

si s n (9(£,£))>o. 



(72) 



Note that in the limit Q — >• 0, x goes to 1 only outside the ergosphere. This does not contradict 
the proposition 0] as it is formulated under the assumption g(d/dt 7 d/dt) < 0. In the region 
g(d/dt,d/dt) > 0, x tends to when Q approaches zero. A numerical study of q<jKx for some 
selected cases is shown in figure [5] 
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Figure 4: In these figures we show a representation of the quality factor q = q = kx for the 
Kerr-de Sitter solution with parameters M = 1.5, A = 1.32 and Killing vector £ = d/dt. In the 
left picture A = 0.1 and in the right picture A = 1. We have represented the quality factor in the 
plane y — z by contour curves (contour interval 1/1000) and only curves with values greater than 
99/100 are represented (this is the region we regard as "close" to unity). The coordinates have 
been compactified through the relations y — ¥ tany, z — > tanz. In the graphs, y G (— tt/2, 7r/2), 
z G (0, 7r/2) (the part with z G (— tt/2,0) is obtained by a reflection under the axis z = 0). We 
have also included a representation of the set of points in which the Killing vector d/dt is timelike 
(shaded region) and the points where the determinant spanned by the Killing vectors d/dt, d/d<j) 
is vanishing (vertical lines and horizontal line). There is a region where the quality factor is above 
999/1000 which is bounded by the thick contour line (it contains the singularity y = 0, z = as 
an accumulation point). This region is larger in the left picture as A is smaller in this case. One 
can check that the size of this region increases when A approaches zero. 



6.5 The Lense-Thirring approximation 

The Lense-Thirring metric can be regarded as a first order approximation to the (exterior) metric 
of a rotating compact body. The metric tensor components arc obtained from (|2ip by just 
retaining the leading terms 

- l^ 2M - -A (7^ 

9tt — — J- H — , 9tx i — ~ e ijk — -Jj, — > 9x i x' ~ 

Without loss of generality we may adopt the choice S 1 = S 2 = (i.e. we choose a coordinate 
system adapted to the body rotation axis). For our computations, it is convenient to express the 
above metric in spherical coordinates (R,9,4>) which are related to the coordinates (a; 1 , x 2 , x 3 ) 
through the usual relation x 1 = R sin 9 cos 0, x 2 = R sin 8 sin 0, a; 3 = Rcos9. We present in figure 
|5] a numerical study of the quality factor q a q for a certain choice of M and of the non- vanishing 
component S 3 . 



7 Conclusions 

Taking the invariant characterizations of the Kerr solution written in terms of the Mars-Simon 
and the space-time Simon tensors we have constructed a number of dimcnsionlcss, non-negative 
scalar quantities (quality factors) which attain the value +1 in a Ricci flat space-time if and 
only if it is locally isometric to the Kerr solution. These quality factors depend on the existence 
of a time-like Killing vector in the space-time and as such they can only be used on stationary 
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Figure 5: Contour curves for the quality factor qq a = qq c = q a n>c for the Kerr Newman solution 
with £ = d/dt. The coordinate R (horizontal axis) has been compactificd with the transformation 
R — > tan R and the coordinate 9 runs along the vertical axis (hence the upper and lower parts of 
the plots represent the axis of the Killing field d/d(j> and they should be identified). The values 
of the mass and angular momentum per unit mass are M — 0.4, A — 0.35 and we consider two 
different values of the charge parameter Q. On the graph of the left we set Q = 0.19 and in 
this case there are two Killing horizons represented by the two parallel vertical lines. The dashed 
curves are the ergosurfaces and the dot is the ring singularity located at R = 0, 6 = 0. Only 
contour curves with values of the quality factor greater than 99/100 have been represented (the 
contour interval is 1/10000). On the graph of the right Q = 4 and there are neither Killing 
horizons, nor ergosurfaces. but the central dot representing the ring singularity remains. In this 
case the contour interval is 1/100 and no further restrictions on the contour curves have been 
placed. In both cases the quality factor approaches unity for R — > ±ir/2 (the asymptotically flat 
regions). 

space-times. However, we must stress that in all the examples studied in section [6] the positivity 
properties of the quality factors were kept even in those points where the Killing vector is not 
time-like. This strongly suggests that all the quality factors studied in this paper are indeed 
non-negative regardless of the causal character of the Killing vector entering in their definition. 
A rigorous proof of this fact has eluded us so far. 

Another interesting question regarding the quality factors is whether they fulfill some hyper- 
bolic differential equation similar to eq. (fT9|) for the Mars-Simon tensor. If that was the case 
then one could analyse the Cauchy problem of such an equation and the global existence of its 
solutions. This could serve to address the non-linear stability of the Kerr solution under the 
assumption of the existence of a Killing vector. 

We have already mentioned in the introduction the uniqueness results dealing with the Kerr 
solution and the recent progress in this direction. In pQ a uniqueness result for the Kerr solu- 
tion was proven which applies to smooth asymptotically flat stationary vacuum solutions of the 
Einstein equations. Besides the usual assumptions of the space-time being the exterior region of 
a regular black hole, a hypothesis about the space-time being "close" to Kerr is assumed. This 
"Perturbation of Kerr assumption" states that the components of the Mars-Simon tensor with 
respect to an orthonormal frame are bounded by a certain quantity when evaluated on a partial 
Cauchy hypcrsurface of the domain of outer communication. One could entertain the possibility 
of expressing this assumption in terms of one of the quality factors introduced in this article. 
This might result in alternative, maybe more powerful, uniqueness results for the Kerr solution. 

Finally, there is the open question of how to generalize the quality factors to the generic 
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Figure 6: The graph of the left is a plot of the contour curves for the quality factor q a q for 
the case of the approximate Lense-Thirring solution with parameters M = 2, S 1 = S 2 = 0, 
S 3 = 14 and fixed 8 = 7r/4 (t is supressed so only the R — <p dependency is kept). As in some 
of the previous plots, the contour interval is 1/100. Note that the curves are all vertical lines 
and one can check that this behaviour is maintained for different values of 8, suggesting that the 
contour hypersurfaces have axial symmetry, as expected on intuitive grounds. We represent one 
of these contour hypersurfaces in the graph on the right. In this graph the coordinates x 2 and 
t have been supressed and the value of q a q is 7/10. Due to the axial symmetry one just has to 
rotate the curve drawn in the graph around the vertical line to obtain the corresponding contour 
hypcrsurface. Note the presence of a region in which numerical value of q a q is less than 1/100 
(the inner boundary of this region is R = 2M = 4). 

situation without Killing vector fields. 
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